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Abstract. A generalisation of the Ginibre ensemble of non-Hermitian random square 



(N 
> 

Qs ' matrices is introduced. The corresponding probability measure is induced by the 

ensemble of rectangular Gaussian matrices via a quadratisation procedure. We derive 

^^ \ the joint probability density of eigenvalues for such induced Ginibre ensemble and 

study various spectral correlation functions for complex and real matrices, and analyse 
universal behaviour in the limit of large dimensions. In this limit the eigenvalues 
of the induced Ginibre ensemble cover uniformly a ring in the complex plane. The 
real induced Ginibre ensemble is shown to be useful to describe statistical properties 
of evolution operators associated with random quantum operations, for which the 

, J dimensions of the input state and the output state do differ. 

rS 

1. Introduction 

In 1965 Ginibre introduced a new three fold family of non-Hermitian Gaussian random 
matrix ensembles as an extension to the mathematical theory of Hermitian random 
matrices [1]. Since then Ginibre ensembles of random matrices have found numerous 
applications in various fields of physics. They can be used to describe non-unitary 
dynamics of open quantum systems [2], transfer matrices in mesoscopic wires [3], 
directed quantum chaos |4]. Real non-symmetric random matrices can be used in 
mathematical finances to describe correlations between time series of prices of various 
stocks [5l[6], and in physiology to characterize correlations between data representing 
the electric activity of brain [71IH]. The same ensemble of real Ginibre matrices describes 
spectral properties of evolution operators representing random quantum operations [9] 
and is useful in telecommunication to characterize scattering of electromagnetic waves 
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on random obstacles [TU], and to describe the effects of synchronization in random 
networks [TT] . 

In his original paper Ginibre derived the eigenvalue distribution of matrices with 
i.i.d. real, complex or quaternion-real entries [1]. These Ginibre ensembles are sometimes 
denoted in literature [I2] as GinOE, GinUE and GinSE, respectively. The letter U, O and 
S stands for orthogonal, unitary and symplectic symmetry class. The case GinOE of real 
asymmetric matrices proved to be the hardest and Ginibre studied only the special case 
that all eigenvalues are real. 

It took another 25 years for Lehmann and Sommers [13] and Edelman [13] to derive 
the complete distribution of eigenvalues for the real Ginibre ensemble. Further difficulty 
arose in the computation of the eigenvalue correlation functions. In 2007 Akemann 
and Kanzieper succeeded in expressing the smooth complex correlation functions as 
Pfaffians [I5], whereas Sinclair presented a method for averaging over the real Ginibre 
ensemble in terms of Pfaffians [12] . Finally, Forrester and Nagao were able to determine 
the real-real as well as the complex- complex correlation functions as Pfaffians using the 
method of skew-orthogonal polynomials [T6|IT7] . while Borodin and Sinclair gave the real- 
complex correlation in addition to a thorough asymptotic analysis [IE] . Simultaneously 
and independently, Sommers [19] and Sommers and Wieczorek [20] derived the complex- 
complex, real-real, and complex-real eigenvalue correlation functions via free-fermion 
diagram expansion. A similar progress was made for the chiral real Gaussian ensemble 
[2T| (see also [221 [23] ioi the chiral complex and quaternion-real Gaussian ensembles) 
and for two non- Gaussian ensembles of real asymmetric matrices [2U|25]. A general 
review on non-Hermitian random matrices can be found in [26] , while a recent overview 
on the Ginibre ensembles is provided in [27] . 

A square matrix A of size A^ pertaining to the complex (real) Ginibre ensemble 
consists of A^^ independent complex (real) random numbers drawn from a normal 
distribution with zero mean and a fixed variance [28]. Normalizing to {Ti AA') = N^, 
the joint probability density function of the matrix entries in this ensemble is 

PGiA)^exp(-^TTA^A] , (1) 

where (3 = 1 for real matrices and (3 = 2 for complex matrices. A"^ is the Hermitian 
conjugate of A which for the real matrices is simply the transpose of A. Note that the 
same formula can be used also to define an analogous probability measure in the space 
of rectangular matrices A. The spectral density of the Ginibre ensemble is described by 
the Girko circular law [29] according to which the eigenvalue distribution for large A^ is, 
in the leading order, uniform in a disk about the origin in the complex plane [301432] . 

The Ginibre ensemble was generalised by Feinberg and Zee [331 [31] who studied 
random complex matrices with joint probability density of matrix entries 

PFz(^)cxexp(-Tr V{A^A)) , (2) 

where ^^(A^A) is a polynomial in A^A. It was found in [33l[M] that in the limit of large 
matrix dimensions the spectral measure of ([2l) can only be supported by a single ring (or 
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disk) on the complex plane. A mathematical proof of the "Single Ring Theorem" was 
recently provided by Guionnet et al. [35] who considered random matrices A = UDV 
with independent Haar unitary U and V, and real diagonal D, see also [36l|37] for an 
alternative approach. The Feinberg-Zee ensemble ([2]) fits into this scheme on making use 
of the singular value decomposition in (jj]). This is because the corresponding Jacobian 
depends only on the singular values of A, see, e.g., [38] • The matrices considered 
in [36|[37] are of the form A = UD with Haar unitary U and diagonal -D, however, 
as far as eigenvalues are concerned the two ensembles UDV and UD are equivalent 
because of the invariance of the Haar measure. Yet another approach to obtaining the 
spectral measure of products of random matrices was proposed in [391H0] and extended 
to weighted sums of unitary matrices in [H] 

In this work we study a generalisation of the Ginibre ensemble oi N x N matrices, 
complex (/3 = 2) or real (/3 = 1), specified by the probability distribution with density 

p!2g(G) = Cf (det G^G) ' " exp (-| Tr G^G^ (3) 

with the normalization constant 

\2! J-.-i r(f + L)) ^' 

Here L is a free parameter such that L > 0. Formally, the random matrix ensemble 
([3]) is a special case of the Feinberg-Zee ensemble ([2]) corresponding to the potential 
V{t) = —^{t—L logt). This specific choice of potential makes the model exactly solvable 
in that one is able to obtain the joint probability density function of eigenvalues and, 
consequently, study the eigenvalue statistics to the same level of detail as in the original 
Ginibre ensemble. For example, the distribution of the real eigenvalues for real matrices 
and all eigenvalue correlation functions become accessible. Obtaining this information 
by the techniques used to study the Feinberg-Zee ensemble seems to be hardly possible. 

Ensemble ([3]) will be called induced Ginibre ensemble as we will show in Section 
|2] that a random matrix G from this ensemble can be generated out of an auxiliary 
rectangular Gaussian matrix X of size (A^ + L) x A^ and a random unitary matrix 
distributed according to the Haar measure on the unitary group. Thus the Gaussian 
measure on the space of (A^ + L) x N rectangular matrices is used to induce another 
measure in the space of square N x N matrices. A similar construction is used to 
generate ensembles of quantum states [12] , as the Haar measure on the group of unitary 
matrices of size kN induces a probability measure in the space of mixed quantum states 
of size A^. Random square matrices of size A^ from the induced ensemble can also be 
obtained by means of quadratisation ([5]) of rectangular (A^ + L) x A^ random matrices. 

The induced Ginibre ensemble of real non-symmetric matrices can be defined in an 
analogous way. We study spectral density and eigenvalue correlation functions for both 
induced ensembles of random matrices. We show that the induced Ginibre ensemble 
exhibits interesting behaviour in the limit of large matrix dimensions: Its eigenvalues 
spreading across an annulus in the complex plane, as opposed to the circular law. 
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Secondly, on the level of correlation functions our asymptotic analysis reveals Ginibre 
correlations, supporting the universality conjecture. 

This article is organized as follows. In Section |2] we introduce a procedure for 
quadratisation of rectangular matrices and we derive the induced Ginibre distribution. 
In Section [3] the induced Ginibre ensemble of complex matrices is investigated, its 
spectral density is derived and eigenvalue correlation functions are analysed. Section H] is 
devoted to the induced Ginibre ensemble of real matrices. In Section[5]it is demonstrated 
that the real induced Ginibre ensemble can be linked to evolution operators associated 
to generic complementary quantum operations, which modify the dimensionality of the 
quantum system. Proofs of technical results are relegated to the appendix. 

2. Quadratisation of rectangular matrices and induced Ginibre ensembles 

Consider complex or real rectangular matrices X with M rows and A^ columns. We 
shall assume that M > A^, so that our rectangular matrices are 'standing'. Let Y and 
Z denote, respectively, the upper square block of size N x N and the lower rectangular 
block of size (M — N) xN oi the matrix X. Since the standard definition of the spectrum 
does not work for non-square matrices we provide a unitary transformation W G U{M) 
intended to set the lower block Z to zero : 



w^x = W^ 



Y 
Z 



G 





(5) 



One can easily find such transformations. Assuming that the matrix X has rank A^, 
consider the linear span S of the column- vectors of X. Let qi, . . . Qat be an orthonormal 
basis in 5, and Qn+i, ■ ■ ■ Qm be an orthonormal basis in S-^, the orthogonal complement 
of S in C*^. If we set W = [qi ■ ■ ■ Qm] then (|5]) holds. Obviously, all other suitable 
unitary transformations are obtained from this W by multiplying it to the right by 
the block diagonal unitary matrices diag[?7, V] where U and V run through the unitary 
groups U{N) and U{M — N), respectively. Multiplying W by diag[lAr, y] corresponds 
to choosing a different orthonormal basis in iS"*", and multiplying W by diag[t/, 1m-n] 
corresponds to replacing matrix G by UG. 

It is straightforward to check that any unitary matrix W G U{M) can be 
transformed to the block form 



W 



where C is A^ x (M - A^) , (6) 



(1^ - CCt)V2 c 

-Ct {iM^M-G^Gfl^ 

by multiplying it to the right by block diagonal unitary matrix as abovqil. 
Correspondingly, we shall seek the unitary matrix ly in ([S]) in the block form ^. This 



X Representation ^ can also be derived by exponentiating matrices of the form A 



A 

-At 



with 



A being N x {M — N). Such matrices A form the orthogonal component of u(A^) ® u(A/ — N) in u(M); 
hence the matrices W — exp(^) provide a natural set of representatives for U{M)/{U{N) x U{M—N)). 
Expanding the exponential in Taylor series leads to ^ with C = A{A^ A)^^/^ sm{{A^ A)^/^). 
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additional condition makes the decomposition in (jS]) unique for ('standing') rectangular 
matrices of full rank. Note that if the matrix X is real then G, W, and, correspondingly, 
C in ([6]) can all be chosen real as well. 

Having settled on the choice of W, we can solve the equation in ([5]) for G and G. 

Lemma 2.1. Let M > N. Suppose that Y is N x N and Z is {M - N) x N, and Y is 
invertible. Then there is a unique M x M unitary matrix W of the form Ij^ such that 
1^ holds. The square matrix G in ^ is given by 

G=h^^^,Z^Z-\ Y. (7) 



A direct calculation proving this lemma is provided in Appendix A 



We now have a procedure for quadratising 'standing' rectangular matrices. Of 
course, in the opposite case of 'lying' rectangular matrices (M < A^) one may apply the 
same procedure to quadratise the transposed matrix X^. Thus, any rectangular matrix 
X can be quadratised by a unitary transformation on its columns (or rows, if the number 
of columns is greater than the number of rows), giving rise to a square matrix G. As G 
is a unique solution of equation (j5]), its spectrum characterises algebraic properties of 
the rectangular matrix X. 

Motivated by studies of quantum operations, see Section [5l we want to explore the 
concept of quadratisation in the context of random rectangular matrices. To this end, 
we will consider Gaussian random matrices X, real or complex, with M rows and A^ 
columns, M > A^, so that the probability distribution of X is specified by the measure 



^v{X) ocexp(--TrA:tA:j |dX| 



(8) 



where /3 = 1 or /3 = 2 depending on whether the matrix X is real or complex, and | dX| 
is the flat (Lebesgue) measure on the corresponding matrix space. 

Note that the result of quadratisation depends on the particular choice of elements 
of VT^X which are set to zero by action of the unitary W in the ansatz ^. For instance, 
one could define another matrix G' by assuming that unitary W brings to zero the upper 
rectangular part of the matrix W^X. However, we are going to apply ([7]) for a random 
matrix X, drawn from a unitary invariant Gaussian ensemble. Thus the statistical 
properties of G and G' will be the same. Therefore for a rectangular random Gaussian 
matrix X one may associate by ([7]) a square random matrix G, which we shall call the 
quadratisation of X. 

The probability distribution for square matrices G ([7]) induced by the normal 
distribution for Y and Z can be obtained directly from ([5]). However, it is useful to look 
at this problem from the Singular Value Decomposition (SVD) perspective, especially 
as one can utilise the known Jacobian of the corresponding coordinate transformation. 

Ignoring a set of zero probability measure, the N x N matrix X'^X has N distinct 
eigenvalues Sj, < si < $2 < ■ ■ ■ < s^, and the SVD asserts that X can be factorised 
as follows 

X = Q T.^'^P\ (9) 
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where S = diag(si, . . . , s^), and Q and P are, respectively, M x N and N x N matrices 
with orthonormal columns, so that Q'^Q = P'^ P = In- The columns of P are in 
fact eigenvectors of X^X and, hence, are defined up to phase factor (or sign for real 
matrices). To make the choice of P unique, we shall impose the condition that the 
first non-zero entry in each column of P is positive. Consequently, the matrix Q is also 
defined uniquely via Q = XPS~^/^. Thus, factorisation ([9]) introduces a new coordinate 
system [Q, P, S) in the space of rectangular matrices. In the new coordinates [38] 

du{X) = di2{Q)djl{P)da{^) , (10) 

where dfi{Q) is the normalised invariant (Haar) measure on the (Stiefel) manifold of 
complex (or, correspondingly, real) M x N matrices with orthonormal columns, dfl{P) 
is the normalized measure defined by the maximum degree form 

(P) = / Ai<k<j<NiMP^dP),k A Im(Ptrfp)^.,) (/3 = 2), 
"^^ ' \ Ai<k<,<N{P^dP),k (/3 = 1), 

on the manifold of unitary (real orthogonal for 13 = 1) N x N matrices satisfying the 
column condition above (the first non-zero entry in each column is positive), and 

N 

da(S) oc (detS)i(^^^^+i^|) ^-^xrE J] |,, - s,f Hds, • (11) 

j<k j=i 

It is apparent from (fTOj) that the matrices Q, P and S are mutually independent. 

Let us introduce an additional unitary (real orthogonal for /3 = 1) matrix U of size 
N X N and rewrite (Q in the form 

X = QUU^T}''^P^ = QUG. (12) 

The matrix G = WJ^^'^P'^ is N x N. Now choose U to be Haar unitary (real orthogonal 
for /3 = 1) and independent of Q, P and S. Then, by rolling back from fITT]) to (^ with 
Q replaced by W , the square matrix G is distributed according to the measure 

d/iindG(G') OC (detGtG)2(*'^-^)exp ( -^TrG^G j | dG\ , (13) 

with the determinant on the right in flT3l) originating from the one in f fTTj) . This is the 
induced Ginibre distribution introduced in Section [1] Because of the invariance of the 
Haar distribution, the unitary matrix U in ( 1T2|) can be absorbed into Q . In other words, 
we have decomposed the rectangular Gaussian X into the product 

X = QG (14) 

of two independent random matrices: the rectangular matrix Q := QU with orthonormal 
columns, this one has uniform distribution, and a square matrix G, this one has induced 
Ginibre distribution. Decomposition ( 1T4l) can also be written as 

X = W ^^ , (15) 
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where W is an M x M unitary (real orthogonal for /3 = 1) matrix obtained from the 
M X N matrix Q by appending suitable column- vectors. This is nothing else as equation 
([5]). One can transform the matrix W to the block form of ([6]), so that ( 1T5l) becomes 



X 



:iyv-cct)i/2 c 



G 





where G = UG for some unitary U. Obviously, G has the same distribution as G. Thus: 

Lemma 2.2. Suppose that X is M x N Gaussian, M > N , with distribution ^. Then 
its quadratisation G has the induced Ginibre distribution U^] . 

This result holds in the cases of complex and real Ginibre matrices. It shows that 
the notion of quadratisation of a rectangular matrix is specially justified for Gaussian 
rectangular matrices. In this case the statistical properties of the outcome do not depend 
on the particular way, how the quadratisation is obtained out of the initially rectangular 
random matrix. This is in a clear analogy to the ensembles of truncations of unitary [13] 
and orthogonal [25] random matrices, the statistical properties of which do not depend 
on the choice of the rows and columns to be truncated. Thus for any random rectangular 
matrix X we may introduce its quasi spectrum as the spectrum of its quadratisation. 

Lemma 12.21 together with Lemma 12.11 provides a recipe for generating matrices 
from the induced Ginibre distribution starting with Gaussian matrices. Interestingly, 
by rearranging f lT2|) one obtains another recipe which is might be less efficient 
computationally but still interesting from the theoretical point of view. Indeed, since 
(XtX)i/2 = PT}l^P\ it follows from ^ that 

G = U^Q^X = U^P\X^Xy/^ = U\X^Xy/^ , where U = PU . (16) 

Recalling that the Haar measure is invariant with respect to right (and left) 
multiplication, one arrives at the following recipe for generating matrices from the 
induced Ginibre distribution. 

Lemma 2.3. Suppose that U is N x N Haar unitary (real orthogonal for P = 1) and 
X is M X N Gaussian with distribution ^ and independent of U. Then the N x N 
matrix G = U{X'^XY^'^ has the induced Ginibre distribution fT^) . 

Obviously, our arguments extend to random rectangular matrices with invariant 
distributions other than Gaussian, e.g. the Feinberg-Zee distribution with density 

PFz(^)cxexp[-Trl^(XtX)], (17) 

where X is MxN,M>N. On applying the procedure of quadratisation to such an 
ensemble, one obtains the induced Feinberg-Zee distribution 

PindFz(G') oc (detGtG)2(^^-^)exp[-Try(GtG')] . 

Another example is provided by rectangular truncations of random unitary or orthogonal 
matrice4§|. By applying quadratisation, one can extend the study of square truncations 
of Haar unitary [43] and orthogonal [25] matrices. 



§ For a range of parameter values the matrix distribution of such truncations is given by (|17|) with 
V{X''X) replaced by a sum of powers of logX^^X and log(l — X^X), see [44] 
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In the context of eigenvalue maps, it is instructive to embed rectangular M x N 
matrices X into the space of M x M matrices by augmenting X with M — N zero 
column- vectors and write the quadratisation rule (j5]) in terms of square matrices albeit 
with zero blocks: 



W^X = G, with X 



Y 
Z 



and G 



G 




where as before W is M x M unitary, Y and G are N x N and Z is (M — N) x N. 

By construction, zero is an eigenvalue oi X of multiplicity M — N and the remaining 
N eigenvalues of X are exactly those of its top left block Y. Thus our quadratisation 
procedure induces an eigenvalue map: the eigenvalues of X are mapped onto the 
eigenvalues of G. Under this map, the zero eigenvalue stays put, and its multiphcity is 
conserved, and, otherwise, the eigenvalues of Y are mapped onto those of G. 

For Gaussian matrices the eigenvalues of Y, for large matrix dimensions (M ^ 1 
and A^ oc M) are distributed uniformly in a disk, and we shall show in the subsequent 
sections that the eigenvalues of G are distributed uniformly in a ring. This hole in the 
spectrum created by quadratisation can be interpreted as due to repulsion of eigenvalues 
of G from its zero eigenvalues. 

3. Complex induced Ginibre ensemble 

The complex induced Ginibre ensemble is defined by the probability density PiaddG) 
(E]) - dll) on the set of complex N x N matrices G = [gjk] with volume element 
IdGI = n -^^j^ dRefyfjfcdlmgfjfc. In the subsequent analysis we restrict ourselves to 
non-negative integer L, in which case the obtained results can be interpreted in the 
context of quadratisation of M x A^ matrices, with L = M — N being a measure of 
rectangularity as was discussed in Section |2l However, our analysis extends almost 
verbatim to real non-negative L, see |15]. 

The symmetrised joint probability density function P(Ai, . . . , X^) of the eigenvalues 
in the complex induced Ginibre ensemble is obtained from that in the Ginibre ensemble 
p!|[28] by multiplying through by det(GG"'')'^ = Ili^i l-^jP^ ^^^ re-evaluating the 
normalization constant. This yields 

N N / ^ \ 

i=i ^-^ ^ j<k j=l \ j=l / 

Consequently, the n-eigenvalue correlation functions 

i?„(Ai, . . . , A„) = — r- / P(Ai, . . . , XN)d Xn+i ■ ■ ■ d A/ 

{N -ny. J 

follow in the determinantal form, 

i?„(Ai,...,AO = det(irjv(Afc,A0)" , (20) 

V / k,l=l 



\N- 
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via the method of orthogonal polynomials in the same way as for the complex Ginibre 
ensemble [28]. The kernel Kf^[\k, A;) is given by 

/..(A.A,)4.-*-"-™'Ef[^. (21) 

j=o ^•' ^ ^ > 

It should be noted that the right-hand side in (12T!) can be expressed as a difference of 
two incomplete gamma functions. This makes the asymptotic analysis of the eigenvalue 
statistics for large matrix dimensions straightforward. 

Of special interest is the mean eigenvalue density PnW '■= -Ri(A) = Kn{\, A), 

where 7(0, z) is the lower incomplete gamma function |46j . 

j{a,z) := / f-^e-^dt . (23) 



JO 

In the limit when A^ is large and L = Na, a > (which corresponds to quadratisation 
of 'standing' rectangular matrices of size (A^ + L) x N) the eigenvalue distribution (in 
the leading order) is uniform and supported by a ring about the origin with the inner 
and outer radii rj„ = a/L and rout = VL + A^, respectively. More precisely, 

lim pn{VNz) = - \Q{VaTT ~ \z\) - 0(v^ - \z\)] , (24) 

where is the Heaviside function, 0(x) = 1 if x > 0, 6(x) = if x < and 6(0) = |. 
Close to the circular edges of the eigenvalue support, for every angle cj), 

lim p^iirout + Oe"^) = hm p^((r,„ - ^e'"^) = ^ erfc(y20 , (25) 

where erfc(x) is the complementary error function, 

2 r°° 

erfc(x) := -^ / e-''dt . (26) 

V^ Jx 

Hence the eigenvalue density falls from I/tt to zero at a Gaussian rate at the inner and 

outer boundaries of the eigenvalue support. As was observed in [37] the scahng law (125|1 

is universal, see Eqs. (42)-(43) in [37] . 

The scaling limits of the eigenvalue correlation functions -R„(Ai, . . . , A^) in the 

induced Ginibre ensemble can also be obtained from fl20|) - fl2T]) . Not surprisingly in 

the regime when A^ — >■ 00 and L oc N one recovers the same expressions as in the 

Ginibre ensemble [HIIHIIIT] , both in the bulk and at the circular edges of the eigenvalue 

distribution, see [Appendix C 



A quick inspection of the joint probability density function (jpdf) of the eigenvalues 
( 1T9|) convinces that the determinantal power in front of the Gaussian weight in ([3]) repels 
eigenvalues from the origin. When the rectangularity index L grows proportionally with 
A^, the eigenvalues are actually displaced from the origin, resulting in the creation of a 
ring of eigenvalues, see Fig. [H 
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1 -0.5 0.5 1 re(z) 



Figure 1. Spectra of matrices pertaining to induced Ginibre ensemble of complex 
matrices for dimension N = 128 and a) L = and b) L = 32. Each plot consists of data 
from 128 independent realisations. The spectra are rescaled by a factor of 1/y/N + L 
to be localized inside the unit disk. The circles of radius: Hn = \/L/{N + L) (inner 
one) and rout = 1 (outer one) are depicted to guide the eye. 



In the context of the augmented matrix G flTSj) the creation of the hole in the 
spectrum can be interpreted as due to repulsion from zero eigenvalues. There are M — N 
eigenvalues of G located at the origin and the remaining eigenvalues are exactly as those 
of G. When L increases, the zero eigenvalues become dominant and repel the rest of the 
eigenvalues strongly, which leads to the creation of the hole in the spectrum. Looking 
at the polynomial on the right-hand side in f l22|) . it starts with the power |Ap^ and ends 
with the power \\\'^(^+^^'^) = \\\'^^^'^~^\ Also in the 2— dim plot in the complex plane 
we see that (M — A^) eigenvalues have been taken off from the centre near the origin 
and put in top of the surface of the circle of support. Thus the outer radius of the 
support of the density increases from y/N to vM while there is an inner hole of radius 
y/I = VM-N - see Fig. [H 

Now we shall explore a different regime when the rect angularity index L = M—N <^ 
A^. This corresponds to the quadratisation of almost square matrices. In the vicinity of 
the origin the corresponding large- A^ (or equivalently large-M) limit can be performed 
by simply extending the summation in fl2Tl) to infinity. For the mean eigenvalue density, 
Pn{X) = Ri{X) it givefjjjl 



lim p7v(A) 



le-|AP 

IT 



E 



r(j + i + 1; 



7r 



r(L) 



and more generally 



N 



lim -R„(Ai 



An) — det \Korigin\^j -I ^k)) j k=i 



with 



K. 



origin\^^j 



(Aj, Afc) 



l7(L,A,Afc) 



TT 



r(L) 



TT r(L) 



t^-^e-'dt 



(27) 



(28) 



II The complex induced Ginibre ensemble is a special case of the more general non-Hermitian randoin 
matrix ensemble studied by Akemann in j48] . Correspondingly, our ((27)) ~p8 )) ) is a limiting case of Eq. 
(4.5) in 135]. Taking the corresponding limit is straightforward for small values of L. 
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At the origin the eigenvalue density vanishes algebraically, Pat (A) ~ ^ r(L+i) &s A — !■ 0, 
uniformly in A^. Away from the origin, the density reaches its asymptotic value l/vr 
very quickljufl. This plateau extends to a full circle of radius a/ZV: 

lim pn{VNz) = -Q{1~\z\) , 

as in the Ginibre ensemble, and, moreover, for reference points y/Nu, \u\ < 1, one also 
recovers the Ginibre correlations. 

Another quantity of interest is the so-called hole probability A{s) at the origin 
giving the probability that no eigenvalues lies inside the disk Dg = {z : \z\ < s}. For 
finite A^ the hole probability A{s) can be derived from the expression: 

N 



A{s)= /"p(Ai,...,A^)n(l 



XDA^j))d'Xi...d'X 



N 



where XDs denotes the indicator function of D^ by employing the method of orthogonal 
polynomials to yield 



N 



with r(a, x) := / e~'^t°'~^dt denoting the upper incomplete gamma function. In 
the asymptotic regime of almost square matrices taking the large A^ limit, while 
keeping L fixed, results in the easily accessible expression for the hole probability 

4. Real induced Ginibre ensemble 



The real induced Ginibre ensemble is defined by the probability density PindG(^) 

N 



Sj) on the set of real N x N matrices G = [gjk] with volume element I^GI = Yl- ^^-^ dgjk- 



In the following we restrict ourselves to even dimension A^. 

4.I. The joint distribution of eigenvalues 

The difficulty in deriving the joint probability density function for real asymmetric 
matrices is due to the fact that there is a non-zero probability p^^k for the matrix G 
to have k real eigenvalues. In the following it is assumed that G has k real ordered 
eigenvalues: Ai > . . . > A^, while / = ^^^ denotes the number of complex conjugate 
eigenvalue pairs xi±iyi, . . . , xi±iyi ordered by their real part. We adopt the convention 
that i/j > for all j . In the case of two complex eigenvalues with identical real part 
the eigenvalue pairs are ordered by the imaginary part. 

As a consequence the eigenvalue jpdf decomposes into a sum of probability densities 
PN,k,i{Xi, ■ ■ ■ , Xk,xi + iyi, . . . ,xi + iyi), corresponding to having k real eigenvalues and 

% A similar behaviour is found in the chiral Ginibre ensemble, see |49] for a discussion in the context 
of gap probabilities. 



Induced Ginibre ensembles of random matrices 



12 



/ pairs of complex conjugate eigenvalues. In order for PN,k,i to be non-zero k must be 
even, as it is assumed that A^ is even. 

The derivation of the eigenvalue jpdf goes as follows [H], see [ISlEOlET] for 
alternative derivations. In order to change variables from the entries of G to the 
eigenvalues of G and some auxiliary variables the real Schur decomposition is employed: 
G = QTQ^, where Q e M.^^^ is an orthogonal matrix, whose first row is chosen to be 
non-negative and the matrix T G M^^^ is block triangular of the form: 

/ -^1 ■ ■ ■ tik ti^kj_i ■ ■ ■ ti^Tv \ 



T 












tk,k+l 



tk,N 
tk+l,N 




V ■■■ Zi / 

Here A is triangular containing the real eigenvalues Ai 
Z is block triangular containing the 2x2 blocks: 



Afc of G on its diagonal and 






6, 



bjCj > 0, bj > Cj and i/j 



\/biC 



3"^] 



on its block diagonal. 

The Jacobian of this change of variable was already computed in 



1^1 



M{^j}j=i,...,k U {xj ± iyj}j=i,...j)\ Ylih - d 



i>k 



with A({2p}p=i^ „) := ni<7('^i ~ ^i) denoting the Vandermonde determinant. 
Consequently we arrive at the relation: 

k I 



P^,GiG)\dG\ = Gl\J\ n |A,I^ n (^- + ^mCn 



m=l 



1 Y^fe \2_i Y^ 

2 ^i=l ^j 2 ^t,- 



^U (- 



2 

2 + i^ 
J^ 2 



X 



dO||dA||dr^||dZ|. 



Here | d O | denotes the volume form on the space of orthogonal N x N 

dT^I = 
In addition 



matrices with positive first row while |dA| = rifcidAj]^ 



N 1 , 

m<n<k'^''mn, 



N 
m=l,...,k,n=k+l,...,N 



dtmn and \dZ\ = U]=i(^bj dcjdxjlll,n=k+i^tr. 



n 

another change of variable is necessary from the entries Xj, bj, Cj of the matrix blocks Zj 
to the real and imaginary part Xj,yj of the complex conjugate eigenvalue pairs and an 
auxiliary variable 6j. The change of variable is performed in the following way: 



Set 



6, 



I (^, + V?+4^ and c, = i [-6, + ^6] + Ay]) , 



(30) 



which implies yj = ^/bjCj and 6j = bj — Cj. The Jacobian of this second change of 
variables can easily be determined: 

4%- 



1^1 



^1 



%| 
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Integrating out the auxiliary variables 6j for j = 1 , . . . , m and t jj as well as using 
Vol(0[A^]) = ^jTv — fTTT finally yields the partial eigenvalue joint probability density 
function: 

^2l-^N(N+l)^-NL 

PN,k,ii^i^ . . . , Afc, xi + iyi, ...,xi + iyi) = ^ ^fi+j-. ^({^i}i=i U {xj ± %}5=i) 

k I 

n |A,l^e-il^^l' n (^™ + yifey--^hmeMV2yJ 

j=l m=l 

where erfc(x) is the complementary error function (1261) . Xj G M for j = 1, . . . ,k and 
Xm + Wm £ C+ for m = 1, . . . ,/. Integrating the partial eigenvalue jpdf PN,k,i over 
R^ X C^' gives pN,k- 

4-2. The {K' , L') — correlation functions 

Again we are interested in the correlations between the eigenvalues. The first starting 
point are the {K', L', k', Z')-partial correlation functions which are just the symmetrised 
marginal probability density functions of K' real eigenvalues and L' complex eigenvalue 
pairs in the case that the number of real eigenvalues is k' while the number of complex 
eigenvalues is /' with different normalisation: 

PN,k',vi.^i^ • • • ' ^k', xi + iyi,...,xv + iyi')dXk'-K'+i ■ ■ ■ dXk'dxv_L'+idyv-L'+i ■ ■ ■ dxvdyv . 

The {K\ L')-correlation functions which are the symmetrised marginals of K' real 
eigenvalues and V complex eigenvalue pairs with different normalisation constant then 
decompose into a disjoint sum of probability density corresponding to having k' real 
and /' complex eigenvalues. They are defined as follows: 

Rk',L'{Xi., . . . ^Xk'^Zi, . . . ,Zli) = 2_^ R{K',L',k',l'){Xl, ■ ■ ■ , Xk', Zi, . . . , Zl') . (31) 

(K'M) 
K'<k',L'<l' 

Remarkably it is possible to express the {K', L')-correlation functions in closed form 
using Pfaffians [T6l[T8l[20l[27]. An elegant approach to the derivation of the correlation 
functions is the method of skew-orthogonal polynomials [iGlfTS] . 



A family {(lj}j=i,..., of skew-orthogonal polynomials is said to be skew-orthogonal with 
respect to the skew-symmetric inner product (—,—), if it satisfies 

iq2j,q2k) = (g2i+i,g2fc+i) = (32) 

(g2j, q2k+i) = -{q2j+i, q2k) = rj6jk for j,k = 0,1,... . (33) 

In the context of the induced Ginibre ensemble (— , — ) denotes the skew-symmetric inner 
product: 

{f,9) := (/,^)k + (/,^)c 
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/OO /'OO 

/ e-5(^'+s'') sgn{y - x)\xy\^f{x)g{y)dxdy 
OO J —OO 

{f.g)c:='2i [ ey'-^\Tic{V2y){x' + yY 

[fix + iy)g{x - iy) - g{x + iy)f{x - iy)] dxdy. 

The method of skew-orthogonal polynomials leads to the following representation of the 
(i^'.L')-correlation functions in terms of Pfaffians. A detailed derivation can be found 
in [18]. Denote 



q{w) := e a"' w^ y eTic{\/2lm.{w))q{w) 
1 



Tj{w) :-- 



sgn{y-w)qj{y)dy, 



iqj{w) sgn(Im(u;)), 



if u; eM 
if u; G C\R 



Then 

Rk',L'{^1, 



Xk', zi, . . . , zl' 



Pfaff 



KjM ( Xj ,Xj') Kn ( Xj , Zm' ) 
KN{Zm,Xji) Ki\f{Zm, Zm'] 



with the 2x2 matrix kernel 



Km(w,w') :-- 



—S]sf{w,w') IS]s[{w,w') + e{w,w') 



(34) 



(35) 



where 



^ 1 
DSNiw,w') = 2 ^ — [q2j{w)q2j+i{w') - q2j+i{w)q2j{w')] 

3=0 ^^ 

2 1 

Sn{w, w') = 2 ^ — [q2j{w)T2j+l{w') - q2j+lT2j{w')] 

3=0 ^^ 

2 I 

ISn{w,w') = 2 ^ — [t2j{w)t2j+i{w') - T2j+l{w)T2j{w')] 



e(w, w') 



3=0 ^ 



- sgn(w — w 



0, 



if w,w' G 
else. 



The indices j and j' in ( l34l) run from 1 to if' whilst m and ??7,' run from 1 to L', so that 
the matrix inside the Pfaffian has the block structure with the top left and right bottom 
blocks being of size 2K' x 2K' and 2L' x 2L', respectively. 

The entries of the Pfaffian kernel depend on the family of polynomials qj which are 
skew-orthogonal with respect to the inner product (— , — ). The direct computation of 
these polynomials is a tremendous task. As a result a different approach is employed in 
order to determine the required skew-orthogonal polynomials and thus the kernel entries 
in (ESD. 
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As already observed in [T^ HH] for the real Ginibre ensemble the following 
relationship for the entry DS of the Pfaffian kernel in fl35|l holds true: 



^" 1 ,1 

2 J^ — [q2j{w)q2j+i{w') - q2j+i{w)q2j{w')] = — 



J=0 ^ 



{w - w') (det {G - wl) det (G - w'I))n-2, (36) 

where (. . .)n-2 denotes the average over the induced Ginibre ensemble of square matrices 
G of size N — 2, and r^ is the normalisation of the A^— th skew-orthogonal polynomial 
as defined in (132|) . 

On integrating out the 'angular' part of G in (1361) . one is left with the integral over 
the 'radial' part of G, see e.g. [50] . 



{det{G - w I) det {G-w' I)) 



N-2 



j=0 



{e,{GG^))N ^2, 
(■) 



[WW 



i\N-2-j 



Here ej{GG^) denotes the j— th elementary symmetric polynomial in the eigenvalues of 
GG^ . The above relation can also be obtained by expanding the product of determinants 
on the left hand side in Schur polynomials in the eigenvalues of G and then averaging 
over the orthogonal group (which is the angular part of G), see [5l] for a similar integral 
over the unitary group. 

The average {ej{GG^))]si-2 can be reduced to a Selberg-Aomoto integral [28]. As a 
result 



4f-l 



1 1 

2 ^ — [q2'j{w)q2j+i{w') - q2j+i{w)q2j{w')] = — 






[W - W) y= 2^ 



rN 



[WW 



'271 



j=0 



m+j + i] 



(37) 



As already observed in [52] the skew-orthogonal polynomials can now be just "read off' 
using the fact that each qj is monic and of degree j by for example differentiating: 



q2j[w} 



Q2j+l[W) 



Q2J+1 



^(2j-M)!9u2i+i 

1 d^^ 



\2j)\du^i 



1 '' 

— iu - w)T{L + N - l)y 



1 ''■ 

^N ^ 



[WU] 



j=0 



TiL + j + l) 



(WU) 



u=0 



j=0 



T{L + j + l] 



u=0 



Hence for j = 1,2,... the following polynomials were found to be skew-orthogonal with 
respect to the skew-inner product (—,—): 

q^^{w) = w^' q2j+i{w) = w^^+' - (2j + L)w^^~' 

In addition to that the first two skew-orthogonal polynomials are given by: qo{w) = 1 
and qi{w) = w. Similarly the normalisation constant can be found by comparison: 

(g2,,g2i+i) = 2v^r(L + 2j + l). 
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Thus the entries of the Pfaffian kernel can be exphcitly determined. A detailed derivation 
of the necessary computation can be found in Appendix B Here we just state the final 
result in Theorem 14.11 below. 
Let 



t{x, z) 



e -2 



erfc(v2Ini(2;)) 



^ 7AW2*-'^^^^^'^'' 



^27r ' ^ ' r(L) ' 

where r(a, z) is the upper incomplete Gamma function, 



(3^ 



and 



r(a, z) 



sn{z,w) 



rN[x,z 



f'-'e-t dt 



N-2 , 



wz) 



J+L 



j=0 

tp[z) sgn[x)2 2 ~^ 2 2z 



'2tx 



J + 1) 

TV+L-l H 2 ^ 2 



1 1^2 



2' 2 



X 



V^N + L-\] 



(39) 
(40) 



Note that sj^i^z^vS) is symmetric in its variables and t(x, 2) and r/vr(x,2;) are not. 

Theorem 4.1. The entries of the complex/ complex (2 x 2) matrix kernel Ki^{z,w) in 
(P^-( r?3]) are given by: 

DSn{z, w) = {w — z)sn{z, w); Sn{z, w) = i{w — z)sn{z, w); ISn{z, w) = {z — w)sn{z, w). 

The entries of the real/complex and complex/real matrix kernels Kj^j^x, z) and Kj^j^z, x) 
in [34^-f3R) are given by: 

DSn{x, z) = {z — x)sn{x, z); DSn{z, x) = — DS'Ar(x, z); 

Sn{x, z) = i{z - x)sn{x, z)] Sn{z, x) = SAr(x, z) + rAr(x, z) + t{x, z); 

ISn{x, z) = —ispf{x, z) — irpf{x, z) — it{x, z); IS^iz, x) = —ISpf{x, z) . 

And finally, the entries of the real/real matrix kernel K^^x.y) in ^^3^- ^3B^) are given 
by: 

DSwix, y) = {y- x)sn{x, y); Sn{x, y) = swix, y) + r^iy, x) + t{x, y); 
1 



ISN{x,y) 



nL) + r(L) + r(L) I' 



7(L + iV- 1,2/2) 



^L+7V-lg|y2 



L + N -1 
-sgn(y)2 



.lix-y)^ l{L + N-l,xy) 

T{L + N-1) T{L + N-1) 



e-^*^t^+^-2rft 



^+Af"|^U +2 2' 2^ -' 



T{L + N- 1] 



L 1„,2\ 



-^*^t^+^-idt-2*-i^^^^ 



r(L) 



e^2* rrft 
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The eigenvalue densities for finite matrix dimensions A^ can be read from the {K', L')- 
correlation functions ( l35l) speciahsing to the (0, 1) and (1,0) cases. Indeed, 



p%{z) := i?o,i(-,^) = PmK^{z,z) = Sr,{z,z) {z e C+) 
is the mean density of complex eigenvalues whilst 

p^(x) := Rifl{x,~) = Pfaff i^Ar(x,x) = S'Ar(x, x) (x G M) 
is the mean density of real eigenvalues. Note the normalisation 

2 / p%{z) d^z + / p^(x) dx = N . 



(41) 
(42) 



Theorem 14.11 now yields the finite- A^ complex and real eigenvalue densities in a closed 
form 

pUx + iy) = J-yeiic{V2y)ey "" }_^ 



3=0 



TU + L + l) 



y eiic{V2y)e^y^ 



2M 



7(L,x2 + r) 7(L + A^- l,x2 + r) 



r(L) 



r(L + A^-i) 



and 



P^(^) 



1 _ , ^2 ^2(i+L) 

'271 



E 

i=o 



r(j + L + i) 



7(L,x2) 7(L + iV-l,x2) 



t(x, x) + rAr(x, x) 



72^1 r(L) r(Ar + L-l) 

where t(x,x) and rn(x,x) are the functions defined in (l38l) - ( 139|) . 



(43) 
(44) 

(45) 
(46) 




-1 -0.5 0.5 1 re(z 



1 -0.5 0.5 1 re(z) 



Figure 2. Spectra of random matrices from real induced Ginibre ensemble for TV = 128 
and a) M = TV = 128 (no hole) and b) M = TV -h 32. Each picture consists of 128 
independent realisations and the spectra arc rescaled by a factor 1/Vm as in Fig. [1] 
In analogy to GinOE we recognise a non-trivial fraction of real eigenvalues. 
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4-4- Asymptotic Analysis 

In this section we shall investigate the real induced Ginibre ensemble in the scaling limit 
when the free parameter L grows proportionally with the matrix dimension A^, which, in 
the language of quadratisation of rectangular matrices, corresponds to tall rectangular 
matrices which are neither skinny nor almost-square. 

The real and complex eigenvalue densities p%{z) (jHj) and p^{x) P6|) are already 
in a convenient form for the asymptotic analysis. The well known limit relation for the 
error function B6]. 



lim ViVu/e^^'"' eiic{^/2Nw) = ^= , 

combined with the saddle-point analysis of the integral in (123|) quickly gives the limiting 
(mean) density of complex eigenvalues. In the leading order, the distribution of complex 
eigenvalues turns out to be uniform in an annulus with the inner and outer radii rj„ = vL 
and Tout = yL + N 1 exactly as in the complex induced Ginibre ensemble. 

Similarly, the saddle-point analysis of each of the incomplete Gamma functions 
in ( H6|) yields the limiting density of real eigenvalues. In the leading order, the real 
eigenvalues in the induced Ginibre ensemble populate two symmetric segments of the real 
axis, [rin.Toutl and [—Touti ~fin\^ with constant density. The theorem below summarises 
our findings and Figs. [2] and |3] provide a comparison of the analytic results versus 
numerical simulations. 

Theorem 4.2. Suppose that L = Na with a > 0. Then: 

(a) In the leading order as N ^ oo, the average number of real eigenvalues in the real 

induced Ginibre ensemble is \/ ^{y/L + N — yL) and the density of real eigenvalues 
obeys the following limiting relation: 

lim p^i^x) = ^= [e(|x| - Va) - e(|x| - VoTT) 

(b) The density of complex eigenvalues obeys the following limiting relation: 

lim p%{^z) = - \q{\z\ -Va)- Q{\z\ - VaTT) . 

On setting L = in the above results one recovers the expected number y^2N/7r 
of real eigenvalues in the Ginibre ensemble [53] together with the uniform densities of 
distribution of real and complex eigenvalues [I4l[53] . 

One can examine how quickly the eigenvalue density falls to zero when one moves 
away from the boundary of the eigenvalue support. At the inner and outer circular 
edges away from the real line, one recovers the same eigenvalue density profile as for the 
complex Ginibre ensemble |47] . 

Theorem 4.3. Suppose that L = Na with a > 0. Then, for fixed C, ^'^ and </> 7^ 0, tt.- 
lim p^((v^-Oe^*) = lim p%{{VlTN + Oe"^) = ^erfc(v^O- (47) 
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It follows from f l47p that the density of complex eigenvalues in the real induced 
Ginibre ensemble falls to zero very fast (at a Gaussian rate) away from the boundary 
of the eigenvalue support. This is also true of the density of real eigenvalues: 

Theorem 4.4. Suppose that L = Na with a > 0. Then, for fixed .^ G M, 



lim p^(yZ-0 = lim pfjiy/TTN + O 



eTic(V2n + ^e-«' erfc(-0 
2^2 



7V-s>oo N-^oo V27r 

The proofs of Theorems I4.2[ 14.31 and 14.41 are straightforward but tedious and are 
omitted. 

Another interesting transitional region appears close to the real line. Here the 
density of complex eigenvalues is more sparse: for finite matrix dimensions p^{x+iy) oc y 
for small values of y. One can easily obtain the complex eigenvalue density profile in the 
crossover from zero density on the real axis to the plateau of constant density far away 
from the real axis. For example, at the origin limAr_j.oo P^iiv) = \/ ^v eric{\/2v)e'^'" , and 
more generally 

lim p^{y/Nu + iv) = \ -veTic{V2v)e'^''^ for |n| G (a/o, V« + !)• 

It should be noted that apart from the support of the eigenvalue distribution 
which clearly depends on a, the limiting eigenvalue density profiles in various scaling 
regimes in the induced Ginibre ensemble are independent of a and coincide with those 
for the original Ginibre ensemble. This correspondence also extends to the eigenvalue 



correlation functions. We show in Appendix D that the eigenvalue correlation functions 



in the induced Ginibre ensemble in the bulk and at the edges are given by the expressions 
obtained for the Ginibre ensemble [18], see also [T6 | fT7 t fT9] . 

4-5. Almost square matrices 

Another interesting regime arises when the rectangularity index L is fixed instead of 
growing proportionally with matrix size as discussed toward the end of Section [31 In 
the bulk, i.e. at a distance of order ^/N from the origin one recovers uniform distribution 
of eigenvalues (real and complex) and Ginibre correlations, whereas in the vicinity of 
the origin new eigenvalue statistics arise. The eigenvalue densities can be obtained by 
extending the summation in (l43l) . ( l45l) to infinity. This yields 



lim p%{x + iy) = \ -yeiic{\/2y) 



2..„„.„../7T.A^2.^ 7(^,x^ + r) 



for the density of complex eigenvalues and 



lim pj^{x 



l{L,X^) , ^_l,2^^2f-l V2'2 



TikAx^ 



T{L) T{L) 



for the density of real eigenvalues. 

As in the case of complex matrices the higher order correlation functions at the 
origin are non-universal: 
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0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1 



Figure 3. Radial density of complex eigenvalues of 256 realisations of a) complex and 
b) real induced Ginibre matrices of dimension iV = 128 with L = 32 compared with 
the analytical results (solid lines) obtained from Eq. (|22|) and (|44| . respectively. Linear 
character of the curves between the inner radius, rin — \fLl\jL + TV ss 0.447 and the 
outer radius, rout = 1, (both rescaled by \/\/L + N and represented by vertical lines) 
suggests that the distribution of eigenvalues on the ring is uniform in both cases. Since 
the dimension N is relatively small, the area on which P{z) is linear does not cover 
the entire interval [7'in,rout]- 



(i) The limiting real/real kernel is given by a 2 x 2 matrix 



^ originy 1 ' ) 



'2ti 



-e"5(--')^2^-t(r,r') 



r{L) 



where 



7(L,r 



I2\ 



T{L) 



+ e 2 



r(L) 



+ I -22"^- 



l(r.-r')2 7(-^5^'^') , /^'^ ea*" „l i"V2'2 



r(L) 



e^H^dt 



(ii) The limiting complex/complex kernel is given by a 2 x 2 matrix 



-'^originy^i ^ ) 



^\/eTic{V2lm{z)) erfc(y2Im(2')) 

'2n 



X 



l,_,')e-^i^-^ril^ 



i\z 



z)e-\^^-^'f--^ [z 



J\„-^(z-z')'^ l{L,zz') 
r(L] 

^ '^ r(L) 



(iii) The limiting real/complex kernel is given by a 2 x 2 matrix 

1 



K. 



origin 



(r,z) 



X 



erfc(A/2Im(z)) 



(.-r)e-i(-)^^ 

-h{r-z)'^ l{L,rz) 
T(L) 



i[z — re 2 



-e 2^ 



-2e 2 



{r-z)^ l{L,rz) 
T(L) 



i(r.-f)2 7(L££) 

r{L) 
■ it(r, 2) 



Nevertheless setting the reference points at a distance of yN away from the origin then 



yields the universal Ginibre correlation functions as given in Appendix D 
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5. Quantum operations and spectra of associated evolution operators 

We start this section witli a brief review of quantum maps wliicli act on the set of 
density operators. Later on we identify a class of maps, for which the associated one- 
step evolution operator is represented by a rectangular real matrix. We argue that for a 
generic random map from this class such an operator may be described by the induced 
Ginibre ensemble of real matrices. 

5.1. Quantum maps 

A quantum state acting on a (i-dimensional Hilbert space Tid can be represented by a 
density matrix p of size d. It is a positive Hermitian matrix, p^ = p > 0, normalised by 
the trace condition, trp = 1. Let Aid denote the set of all density matrices of size d. 

Consider a linear quantum map $, which maps the set of density matrices onto 
itself, $ : A4d -> Aid- Such a map is called positive as it transforms a positive operator 
into a positive operator. In quantum theory one uses also a stronger property: a map 
is called completely positive if any extended map, $ 1„ is positive for an arbitrary 
dimension n of the auxiliary subsystem. 

Any completely positive map $ can be written in the following Kraus form, 

k 

p' = ^{p) = Y,ErpE\. (48) 

The Kraus operators {£^i}f=i and their number k can be arbitrary - see e.g. [51] • The 
map ( HHll preserves the trace, trp = trp', if the entire set of k Kraus operators satisfies 
the identity resolution, 

k 

Y,E\E, = 1. (49) 

Any physical transformation of a quantum state p can be described by a completely 
positive, trace preserving map, which is also called a quantum operation. 

Alternatively, any quantum operation can be written as a result of a unitary 
operation acting on an extended system and followed by the partial trace over the 
auxiliary subsystem E, 

p' = $(p) = trE [U{p ® \u){u\)U^] . (50) 

Here a unitary operator U G U{kd) acts on a composite Hilbert space T-Ls ® T-Le, where 
S denotes the principal system of size d, while E denotes an environment of size k. It 
is assumed that the environment is described initially by a pure state, |i^) G He- 

The action of the linear transformation (HH!) can be described by a matrix $ of size 

p' = $p or Prn^ = ^mi^pnu , (51) 

where the summation over repeated indices goes from 1 to d. Since the matrices p 
and p' represent quantum states, so they are positive and normalised, the following 
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matrix with exchanged indices, Dmn := $m/i is also Hermitian, positive and normahsed, 
tiD = d. Thus this matrix of rank k, called dynamical matrix is related to a quantum 
state ? acting on the composed Hilbert space Tid <S) T-Ld- The state <;^, proportional to the 
dynamical matrix, can be expressed by an extended map [51] 

^ = \d = ($®l)|^+)(^+| , (52) 

where \ip^) denotes the maximally entangled state, !'?/;"'") = ^ X]j=i K) ® K)- 

The dynamical matrix D is Hermitian by construction, while the superoperator 
matrix $ is in general not Hermitian. Positivity of the map and the trace preserving 
condition imply, that the spectrum of $ is contained in the unit disk and there exists 
at least one eigenvalue equal to unity. It corresponds to the invariant state of the map. 
The dynamical properties of the map can be described by the modulus of the subleading 
eigenvalue, which determines the behaviour of the system under consecutive actions of 
the map. 

For a generic random map, distributed according to the flat (Euclidean) measure 
in the space of quantum operations, the subleading eigenvalue is almost surely strongly 
smaller than one, which implies that the convergence to equilibrium occurs exponentially 
fast. In fact statistical properties of the evolution matrix $ of a random operation 
can be described by the real Ginibre ensemble [55]. The same is true for operations 
corresponding to quantum djTiamical systems under assumption of classical chaos and 
strong decoherence [9l|56]- The real ensemble is applicable here since the map $ 
preserves Hermiticity and the matrix $ becomes real in the generalised Bloch vector 
representation [51] . 

Observe that independently of the dimension k of the environment the dimension 
d of the output state p' does not change under the action of the transformation (l50l) . 
so that $ : M.d — ^ M.d- However, in several applications one uses a more general class 
of quantum maps, which may change the dimension of the state. Among them one 
distinguishes an important class of complementary maps. A map $ complementary to 
$ can be defined as [571158] 



a = ^{p) = iis[U{p®\v){v\)U^]. (53) 

Note that the only difference with respect to (l50l) is that the partial trace is performed 
with respect to the principal system S. The output state a E M.k describes thus the 
final state of the environment after the interaction with the system described by the 
global unitary matrix U . Thus the complementary map sends the density matrix of size 
d into a state of size fc, $ : Aid -^ M.k- Comparing both maps we see that the role of 
parameters d and k is interchanged: the rank of the dynamical matrix corresponding to 
$ is equal to /c, while for the complementary map $ it is equal to d. 

The complementary map $ can also be described in the form (ISTl) . However, as the 
dimension d of the input state p and the dimension k of the output state a do differ, the 
matrix $c representing the evolution operator $ under the action of a complementary 
map is a rectangular matrix of size k'^xd^. To study spectral properties of such evolution 
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Figure 4. Spectra of operators $ associated with evolution operators for generic 
complementary quantum maps for fixed size d = 14 of the initial system and different 
sizes of the environment, a) fc = 10, b) fc = 14 = rf (which implies no hole in the 
support of the spectrum), and c) fc = 18. The data in each panel are superimposed 
from 8 realisations of random complementary operations. Solid circles represent the 
radii implied by Eq. dSF 



operators one has therefore to go beyond the standard diagonahsation procedure which 
holds for square matrices only. 

5.2. Complementary m,aps and real induced Ginibre ensemble 

We are now in position to investigate spectral properties of evolution operators 
associated with a random complementary map defined by ( l53l) . As such an operator 
is represented by a rectangular matrix, we use formula ([7]) to obtain a square matrix, 
which can be diagonalized in the standard way. 

Random complementary maps were generated numerically by selecting a fixed 
initial state {u) of the k dimensional environment and by plugging into fj53|) a random 
unitary matrix U distributed according to the Haar measure on U{kd). This procedure 
allows one to generate random maps according to the Euclidean measure in the space 
of all quantum operations [55] . 

Fig. m shows the exemplary spectra of $ associated with random complementary 
maps, which transform an initial quantum state of size d = 14 into an output state of 
various dimensionality. In the symmetric case k = d the evolution operator is repre- 
sented by a square matrix, so the spectrum covers the entire disk. In the cases k ^ d the 
evolution operator $ is represented by a rectangular matrix $c- Thus the spectrum of 
its quadratisation $, which corresponds to a matrix from the induced Ginibre ensemble, 
covers asymptotically a ring in the complex plane. The radius of the inner ring depends 
on the difference \k — d\. 



Consider a random complementary map $, which sends the set A^^ of d-dimensional 
states into M.k- The evolution operator can be represented by a rectangular matrix, $c 
of size k"^ X d^. Assume first that k > d. Constraints imposed by the trace preserving 
condition ( H9l) are known to be weak even for relatively small dimensions [9l|55] , so it is 
legitimate to assume that $c can be described by a real random rectangular matrix of the 
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Ginibre ensemble. Thus a quadratisation $ of $c will be described by an appropriately 
rescaled matrix of the induced Ginibre ensemble with M = k"^ and N = d"^. 

To find the scaling factor we need to discuss the normalisation of $c- The squared 
norm, ||$c|P = Tr$J$c is equal to the norm ||-D|p = TrD^ of the Hermitian dynamical 
matrix, which contains the same entries in a different order [51]. Since the map $ 
changes the dimension d of the input state into k, expression ( !52l) implies that the 
corresponding dynamical matrix D is a dk x dk Hermitian matrix of rank d. 

Neglecting the constraints implied by the trace preserving condition, Ti bD = 1, 
we assume that in case of a random complementary map the dynamical matrix D 
behaves as a random mixed state generated by an induced measure [12] . In such a case 
the average purity of a random state p of size A^ obtained from an initially random 
pure state of size NK by the partial trace of an over the K dimensional environment 
reads, (Trp^) = [K + N)/KN. In the case of dynamical matrix of a complementary 
map we need to substitute the correct dimensions, N ^ kd and K ^ d, and use the 
normalisation constant D = dq io obtain an approximate expression 

(TrZ}2) = d\Tie) ^ rf2M±^ = d{k + V)_ ^ ^_ 

d^k k 

This in turn implies the relation (Tr($J$c)) ~ d, so the average norm of the 

superoperator $c representing the complementary channel $ does not depend on the 

parameter k and can be estimated by ||$c|| ~ yd - see Fig. [51 
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Figure 5. Average squared norm of the superoperator related to the complementary 
map, t = (Tr$J$c) as a function of size d of the input system. Various symbols denote 
different values of the size of the output system, fc = 5 (o), fc = 10 (n), k = 15 (x), 
fc = 20 (+) and A; = 25 (•). Dashed line is plotted to guide the eye. It corresponds to 
fc = 25 and shows that the bigger value of fc, the better approximation to the asymptote 
represented by the solid diagonal line. 
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As discussed in previous sections, an M x A^ matrix X (with A^ < M) of the 
Gaussian ensemble (|8]), for which the average norm satisfies {TtX'^X) = NM, leads to a 
square random matrix G of size A^ from the induced Ginibre ensemble ( IT3l) by the way of 
quadratisation. Its spectrum is (asymptotically) supported in the ring r G [rjn, rout] with 
Tin = vL and rout = VL + N. To apply this ensemble for a rectangular superoperator 
matrix <l>c it is then sufficient to substitute M ^ k'^ and N ^ d"^ and to use the rescaling 
^^ = X/kyd to match the normalisation. Therefore the spectrum of the quadratisation 
$ of the superoperator $c associated with the complementary map $ forms a ring of 
the inner radius rjn = a/1 — d^ /k"^ /\fd while the outer radius reads rout = l/^fd- 

Consider now the other case, in which the dimensions determining the 
complementary map satisfy k < d. Then the superoperator $c, represented by a 
rectangular matrix k"^ x d^, gets its quadratisation by the same procedure from the 
transposed matrix as mentioned in Section |2l and again it is described by the induced 
Ginibre matrices with the level density analysed in Section HI The key difference with 
respect to this distribution, describing the transposed rectangular matrix X^, is that 
the role of variables M and A^ is exchanged. Hence in the case for the transposed 
superoperator $^, one needs to exchange the dimensions k and d. Taking this into 
account and the normalisation, ||"^c|P ~ d, one infers that the outer radius of the ring 
reads in this case rout = y<P /kyd = ydjk. 

Thus the square matrices $ associated with the superoperators corresponding to 
random complementary maps $ : M.^ — ?> A^^ can be described with real matrices of the 
induced Ginibre ensemble. Apart of the leading eigenvalue Ai = 1, which is implied by 
the trace preserving condition, the spectra are asymptotically localized in a ring in the 
complex plane. Both radii of the ring read 



in; ' out J 



J_ /i _ d£ J_ 



V^ /l _ fc^ v^ 
fc V d2 ' k 



for k ^ d 

(55) 
for k ^ d . 



The above predictions show good agreement with numerical data obtained for several 
realisations of random complementary maps and shown in Fig. HI To obtain a clearer 
figure the plot is magnified by a factor vd. 

6. Concluding remarks 

Although the specific example of the Ginibre ensemble is completely solved, the theory 
of non-Hermitian random matrices is still far from being as thoroughly understood 
as its Hermitian counterpart. In this work we have introduced a new generalisation of 
the ensemble of non-Hermitian Ginibre matrices and derived explicit results for spectral 
density in the complex and real cases. Using the method of skew-orthogonal polynomials 
we derived various spectral correlation functions which can be expressed as Pfaffians. 
Analysing asymptotic behaviour in the limit of large matrix dimensions we have found 
a universal behaviour of the eigenvalue statistics. 
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The induced Ginibre ensemble of square matrices of size A^ is parametrised by a 
single discrete parameter L > 0. In the case L = the model reduces to the standard 
Ginibre ensemble [1], with the eigenvalues, for large A^, uniformly distributed in the disk 
of radius y/N about the origin, while for L > the eigenvalues are repelled from the 
origin and form a ring in the complex plane vL < |z| < \/N + L. 

This form of the spectrum suggests a comparison with the model ([2]) of Feinberg- 
Zee, for which the 'single ring' theorem was proven [33l[35l|59]. Our model formally 
belongs to the Feinberg-Zee class, with the potential 

V{G^G) = G^G-L log G^G, (56) 

but due to the log function the assumption that the potential is polynomial is not 
satisfied. Analysing the joint probability distribution for induced Ginibre ensemble it 
is easy to see that eigenvalues near the origin are unlikely, which implies the ring of 
eigenvalues. On the other hand it is not simple to find a specific mechanism responsible 
for the single ring distribution for the general version of the Feinberg-Zee model. 
Furthermore, in the latter model the eigenvalue correlation functions are not known, 
while we could compute them for the induced Ginibre ensemble for the complex and 
real versions of the model. 

Our work leads to a straightforward explicit algorithm to generate random matrices 
from the induced Ginibre ensemble. Its is sufficient to take a rectangular M xN random 
Gaussian matrix X, construct the positive Wishart like matrix, X^X of size N, take its 
square root and multiply it by a random unitary U distributed according to the Haar 
measure on the unitary group U{N). The result G = U VX^X is distributed according 
to the desired joint probability function ([3]), while its singular values are described by 
the Marchenko-Pastur distribution with parameter c = M/N. This follows from the 
fact that G^G = X'^X, hence the square random matrix G has the same singular values 
as the initial rectangular Gaussian matrix X. 

Alternatively, square random matrices from the induced Ginibre ensemble can be 
obtained by the quadratisation (I5])-(|6]), ([7]) of rectangular Gaussian matrices. Although 
for a given rectangular matrix various quadratisation algorithms produce different square 
matrices, the statistical properties of square matrices obtained by quadratisation of 
rectangular random Gaussian matrices do not depend on the choice of the algorithm. 

An analogous construction works also in the real case. Take a rectangular M x N 
real Ginibre matrix X, construct the square root of the positive Wishart matrix, X'^X 
and multiplying it by a random orthogonal matrix O distributed according to the Haar 
measure on the orthogonal group 0{N). The result G = OVX'^X is then distributed 
according to the induced ensemble of real Ginibre matrices, specified by ([3]) with /3 = 1. 

The induced ensemble of random Ginibre matrices offers a simple model for further 
research within the developing theory of non-Hermitian random matrices. We are 
tempted to believe it will find its applications in several fields of physics. For instance 
they can be helpful for analysis of evolution operators associated to certain classes of 
generic quantum maps. In the simplest case as the dimension of the input and output 
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states are the same, the associated evolution operator can be described by the standard 
ensemble of real Ginibre matrices [9]. However, studying complementary quantum 
operations, or other quantum maps in which the dimensions of the initial and the final 
states do differ, one copes with evolution operators represented by rectangular matrices. 
Our work provides evidence that the statistical properties of square matrices associated 
with the evolution operators for complementary quantum maps do correspond to the 
induced ensemble of real Ginibre matrices. 
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Appendix A. Proof of Lemma 12.11 

Assuming IV as in (j6]), by multiplying through in ([5]), one obtains an equation for C, 
C^Y + {1m-n - C^CY^^ Z = 0. Hence, 

Z = - (1m_jv - C^Cy^^^C^Y = -C^ {In - CC^y^^^Y . (A.l) 

Consequently, by making use of (^ again, 

G={1n- CC^) '^' Y-CZ={1n- CC^) "'^' Y. (A.2) 

It is easy to check that 

Z^Z = Y^ {In - CCt) "V - Y^Y (A.3) 

which implies that 

yty + z^Z = Y^1n- CC^) '^ Y. (A.4) 

This in turn allows us to write 

(1^ - CC^y' = ^ (yty + z^z) 1 (A.5) 

which when substituted into l\A.2\} yields (jTj). Note that the desired result ([7]) can be 
also rewritten in a more symmetric form, 

G = Y (rty)-i/2 1^^^ ^ (yty)-i/2 ^t^ (Fty)-!/^^/^ ^^ty) V2 ^ ^^ g^ 

which shows that all matrix square roots operate correctly on positive definite objects. 
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In this appendix we determine the finite N correlation functions for the real induced 
Ginibre by determining the matrix entries of the Pfaffian kernel given in [M] and hence 
prove Theorem 14.11 Equipped with the appropriate skew-orthogonal polynomials and 
their normalisation the task of determining the entries of the matrix kernel for the 
(j,m)— correlation functions can now proceed. Firstly (137|) implies for w, 2 G C: 

N-2 



DSm(z,w) 



'2ti 



^(wU(z)(z -w)Y -r^ 7 



(B.l) 



Noting that Sn{z,w) = iDSis[{z,w) and ISj\f{z,w) = —DSn{z,w) we have completely 
determined the entries of the complex-complex matrix kernel. 



Let us next consider the case x G M, 2 G C. The following approach is borrowed 
from [16]. We observe that 

d 



i-r2 -T 
-62^ X ^ 



dx 



e-h\^i+L 



\^\j+L 



q2j+i[x) = 

which implies for j > 0: 

T2j+i{x) = e"5- X- 
Furthermore direct computation shows that: 

Ti{x) / sgn(x - t)e~^^^x^~^dt = e'^^'^x'^^'^^ 

All in all 

N -, AT -I 



tp{z)z^ Ti{x) / sga{x — t)e~2^^x^~^dt 

L 2 Jk 






V27rr(L + 1) 



ih(z)z^- / sgn(x-t)e ^^^x^'^dt 



^T{L + 1) 



1 -iT^-irZ 

-g 2 2 



/27r 



erfc(y2Im(^)) ^ 
i=o 



[xz 



:L + 2j 



1 






r(L + 23—,^V^'^^'^'' r(L + i) 



In addition to that: 



2 1 
S= X] -^2i+i(2;)r2,l 



X) 



j=0 

1 

V2^ 



^(.).-E''^'"^^':"'t''^^'^= + 



i=i 



r(L + 2j + i) 



27rr(L + 1 



-i}{z)z^To{x) 
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Rearranging the summation gives: 

s- 



2ttT{L + N -1 



-tjj{z)z 



L+N-1 



TN-2{x) 



Uw^'E^ 



T2j+2{x) - (L + J + 1)t2j{x)] Z 



2j+l 



21X-' ^ r(L + 2j + i: 



-62^ X 



d_ 
dx 



Another differential equation: 

(l2j+2{x) - (2j + L + l)q2j{x) 
leads to: 

r2M^) - (2J + L + l)r2,(x) = e-'^^'x'^^'^^' 
As a consequence 



^„i,2^2i+L+l 



S 



zil}[z)z 



L+N-lr,^+N 



'Itx 

Finally we obtain 
1 



22- 



2 sgn X 



Af-2 



Sn{z,w) 



e 2^ 



27r 



^wE 



a;^; 



^f:^ + Z^ _ i la;2^ 
' V 2 ^ 2 2' 2 ^ 

r(L + iv-i) 



,L+2i 



^{z)e ^"'E 



i=o 



(^^)L+2j+l 

r(L + 2j + 2) 



i=o 



r(L + 2j + i^ 



1 



^(^)^^+^-i2#+^-isgn(x) ''2 ' 2 



/ \ o ~n Q o 1 o "^ 



2' 2' 



-4{z)T^ 



-1 -*- I 2 ' 2^ J 

^27r'^' - ^ ' r(L + iV-l) ' v^"^~^" r(L + i) 

The last entry requiring explicit computation is ISt^[x^y) for x^y ^^. In this case the 
relationship: 



ISN{x,y) 



SN{t,y)dt 



comes handy. Using the expression obtained for SN{x,y) and in addition to that 
employing the integral representation 

e-^yY: 



(tyY 



j=0 



T{L + 2j + l] 



-f{L,ty) -f{L + N-l,ty) 



T{L) T{L + N-1) 



leads to the following starting point for our derivation 

1 rv 



ISN{x,y) 
+ 



V2ttT{L) 
1 rv 



e 2 



l(^^y)\y-tML,ty)dt 



27rT{L + N - 1] 
^sgn(x)2 2+" 2 



■'^^'-y'^\y-t)-f{L + N-l,ty)dt 



2 2' 2 '' 



27r 

2 2 ^ / 62* t^dt . 



T{L + N-1] 



e^*t^+^-irft 



The above expression can be simplified by employing integration by parts with respect 
to t. As a conclusion we have derived all the possible entries of the Pfafiian matrix 
kernel. 
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In this appendix we state the asymptotic behaviour of the correlation functions of the 
complex induced Ginibre ensemble for different scaling regimes. 

Theorem Appendix C.l (The limiting correlation functions in the bulk). Let 

u,zi, . . . , Zn be complex numbers and set \k = yNu + Zk for k = 1, . . . ,n and L = Na, 
then for u G R: 

~1 



lim Rn(X) = det 



IT 



exp 



Ffc 



+ ZjZk 



j,k=l 



with R= {r e C\y/a < \r\ < Va + 1}. 

Proof. We need to analyze the asymptotics of the kernel: 



KN{Zj,Zk) 



1 -il-/]V«+2j||v^iZ+2 ;l-i 



-e 21 
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,j I — ^ I V Nu-\-Z}^ \'Nu-\-Zj^ 
l+L 



X 



N- 

E 

1=0 



1 (\y/Nu + ZjWVNu + Zk 



T{l + L + 1) 



1 

— I 

71 



;^^-N\u\^-y/N{uZk+uzj)-ZjZk^ — -^ ^ \-ZjZk^^iuZj-uZj)^-^{uZk-uZk) 



^/N / 



Vn t 



N- 



X 



E 



1 {N\uY + VN{uZk + uzj) + ZjZkj 



l+L 



1=0 



r(/ + L + i) 

In order to simplify this expression we resort to the following trick outlined in [18]. We 
set: ^]\[{zj) = e^"*-"^^'""^^'-' and define the diagonal matrix 

Dn = diag NjN{Zl),ijN{z2),- ■ ■ j-ipNiZn) 

We note the following: |'?/'Ar(s)| = 1 and as a consequence 

Rn{\) = det{Dr,)Rn{\)det{Dl). 
Thus we can write for the n-point correlation function 

R (A) = det I" — e^~2 ^-^+ZjZk^-N\u\'^~y/N(uZk+uZj)-ZjZk 

LIT 



2 """^J ^k g 



Af- 



X 



E 

1=0 



I (n\u\^ + y/N{uZk + uZj) + ZjZk) 



l+L 



T{i + L + r 



j,k=i 



Applying another saddle-point analysis it can be shown that for u & R: 

l+L 



N- 



lim ^-^\'^\'^-^i'^^k+uzj)-zjZh V^ 

l/"— inn ^-^ 



N^oo 



1 \N\u^ + VN{uZk + UZj) + ZjZkj 

r(/ + L + i) "^' 



which gives the limiting expression for the ra-point correlation functions. 



n 
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Employing the same technique again, albeit with a slightly different saddle-point 
method, yields the limiting behavior of the n-point correlation functions around the 
circular edges of the eigenvalue density. 

Theorem Appendix C.2 (The limiting correlation functions at the edges). Let 

u,Zi, . . . , Zn be complex numbers with |m| = 1, 

1. Setting Xk = ^/N{a + l)u + z^ for k = 1, . . . ,n leads to the limiting correlation 

functions at the outer edge yjN{a + 1).' 

1 



lim Rn{\i,...,K 



det 



N^oo 



271 



exp 



\Zk\ 



+ ZjZkj ( erfc ''^■^ 



ZjU + ZkU, 



v^ 



hk=l 



The same limiting expression is found around the inner edge yWa of the eigenvalue 
density by setting Xk = yNau — Zk for k = 1, . . . ,n. 



Appendix D. Limiting correlation functions for the real induced Ginibre 
ensemble 



Throughout this section it is assumed that L = Na. 

Theorem Appendix D.l (The limiting correlation functions in the bulk). Let n G M 

such that ^JoL < \u\ < \/a + 1 and let ri, . . . ,rj G M as well as Si, . . . , Sm € C+\]R. 
Furthermore set Xt = yNu + Vt for t = 1, . . . ,j , z^ = yNu + s„ for v = 1, . . . ,m and 
L = Na, then: 



lim Rj^rnixi, ...,Xj,zi,..., Zm) = Pfaff 

TV— >oo 



K{rt,rt>) K{rt,s^i) 
K{s^,rt') K(s„,v) 



where t, t' = 1, . . . , j and v,v' = 1, . . . ,m. 
(i) The limiting real/real kernel is given by 
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(a) The limiting complex/ complex kernel is given by 
1 
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(Hi) The limiting real/complex kernel is given by 
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Proof. In the following only the derivation of the limiting behaviour of the real/real 
kernel will be outlined, as this is the most involved computation. All other results can 
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be deduced in a similar fashion. We start by computing the asymptotic behaviour of 
Sn- Using the results from Appendix C.l it is obvious that: 

lim sn(VNu + r, VNu + r') = ^=e-5('-'^')'. 



The next term to be analyzed is: 

1 2t("+i) 

Tt^iyNu + r, viVu + r') = %<gi\.{yNu + r 



27r ^ ' ^r(iV(a + l)-l) 

We can apply the duplication formula for the gamma function: 



r(2^) = r(.)r(. + V^-i^. 

2 a/tt 



and obtain: 

r^{^fNu + r, ViVn + r') = sgn(yiVn + r)2f ("+i)-te-^""'(l + ^^)A^("+i)-i 



r(y(a + l))r(y(a + l- 2)) V2 2 2 

Furthermore the use of the Stirling formula: 



as well as: 
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leads to: 
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Now it can easily be seen that: 

hm e'^^^''^ " 'e ^ ^ -^ ' = 0. 

W->-oo 

In addition to that for \u\'^ < a + 1 the expression a + 1 — u'^ + 2(a + 1) ln(^^) is 
negative and a saddle point method shows: 

£"L r(f (a + 1-1) " (f'° ^ ^' 4- 5'^" + --'O - "■ 

As a conclusion we have shown that limAr_j.oo rj\f{\^u + r, viVn + r') =0. The next 
term 

tf^iVn + r, TiVn + r') = ^ / , e~^(^"+"'^'(ViVn + r')^°r(— a, -(ViVn + r?) 

^/2^^ T{Na) ^2 2^ ^ ^ 
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can be dealt with in a similar fashion. We have now determined the scaling limits of 
the real/real entries Sn and DSn. The scaling limit for the entry IS^ can be found by 
applying a saddle point method on each of the eight integrals. In addition to that the 
asymptotic relationships derived for rj^j and t can be applied. D 

Similar calculations lead to the conclusion that in the complex bulk and also at the 
edges the eigenvalue correlation functions in the real induced Ginibre are exactly the 
same as those in the real Ginibre ensemble. We omit the derivations and only state the 
results in the two theorems below. 

Theorem Appendix D.2 (The limiting correlation functions in the complex bulk). 

Let u be a complex number such that u & R = {r & 'C\^/a < |r| < \/a'^\^} and let 
Si, . . . ,Sm G C Furthermore set Zj = yNu + Sj for j = 1, . . . ,m and L = Na, then for 
u e R: 



lim i?o,m(-, Zi,...,Zm) = - det 



\Sk\ Sfc/ 
exp ( — + ZkZy 



k,k'=l,...,m 

Theorem Appendix D.3 (The limiting correlation functions at the edges). Let 
u = ±1, ri, . . . ,ri G M as well as si, . . . , Sm G C+. Setting Xj = ^J N{a + l)u + Tj 
for t = 1, . . . , / and z^ = \/N{a + l)u + Sk for k = 1, . . . ,m leads to Ginibre behavior 
for the limiting correlation functions at the outer edge of the eigenvalue distribution as 
described in [7g[ /. In addition at the inner circular edge Xj = yWau — Tj for t = !,...,/ 
and Zk = yNau — s^ for k = 1, . . . ,m the Ginibre limiting correlation functions can 
again be recovered. 
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